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Abstract: In this article, we derive an approximate asymptotic analytical expression for 
the long-time chronoamperometric current response at an inlaid microband (or laminar) 
electrode. The expression is applicable when the length of the microband is much greater 
than the width, so that the diffusion of the electrochemical species can be regarded as 
two-dimensional. We extend the previously known result for the diffusion-limited current 
response (Aoki, K. et al. J. Electroanal. Chem. 1987, 225, 19-32 and Phillips, C.G. J. 
Electroanal. Chem. 1992, 333, 11-32) to accommodate quasi-reversible reactions and 
unequal diffusion coefficients of the oxidant and the reductant. Comparison with numerical 
calculations validates the analytical expression, and we demonstrate that unequal diffusion 
coefficients can substantially change the current response. Finally, we discuss the form of the 
long-time current response for a one-step, one-electron redox reaction if the rate constants 
are modelled in the Butler- Volmer framework, and indicate the importance of choosing the 
width of the microband appropriately to allow accurate experimental determination of the 
standard kinetic rate constant and the electron transfer coefficient. 

Keywords: chronoamperometry; two-dimensional; electrode; band; microband; 
lamina; voltammetry 
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1. Introduction 

Chronoamperometry is a widely used voltammetric technique [1], whereby the potential at an 
electrode is initially stepped from a level at which no current is flowing to a level at which the redox 
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reaction at the electrode can progress and a current starts to flow. The current is generated by the 
exchange of n electrons at the electrode surface according to the following redox reaction with forward 
and backward rate constants, kf and fc& (m-s -1 ): 

Ox + ne^Red (1) 

k b 

The information accessible from the chronoamperometric current response can be enhanced by utilising 
microelectrodes, and in particular ultramicroelectrodes. In addition to reducing the distorting influences 
of double-layer capacitance and ohmic drop on the Faradaic current [2,3], the small dimensions of 
the electrodes allow access to information about fast redox kinetics that would previously have been 
inaccessible [4]. Microelectrodes can be fabricated in many different geometries, such as inlaid disks, 
inlaid ring-disks and mounted hemispheres, but here we consider an inlaid microband [5,6] such that its 
length, L, is much larger than its width, w = 2d, as depicted in Figure 1 . If the length L is allowed 
to be macroscopic, but the width w is kept microscopic, then the overall current response through the 
electrode is increased, whilst the advantages of microelectrodes are retained. From a modelling point 
of view, the advantage of assuming L 3> d is that effects due to the ends of the band can be neglected, 
and diffusion of the oxidant and the reductant can be regarded as two-dimensional (provided that the 
time-scales under consideration are such that t <C L 2 /D, where D is the typical size of the diffusion 
coefficients of the oxidant and the reductant, Do and Dr). By considerations of symmetry, the current 
through an equivalent laminar electrode in free space is twice that of the inlaid band electrode, since it 
must have identical currents flowing through its upper and lower surfaces. 

Figure 1. Cartoon of an inlaid microband electrode of width w = 2d and length L. The 
axial length, L, of the electrode is much greater than the width, w, i.e., L 3> w. Then, for 
time-scales t such that t <C L 2 /D, where D is the typical size of the diffusion coefficients of 
the redox species, this allows the two-dimensional-diffusion approximation to be used, since 
contributions to the current from the three-dimensional diffusion at the ends are negligible. 




4 



w = 2d 



The majority of the theoretical modelling for the current response at an inlaid microband electrode 
has involved numerical simulations, using a variety of different techniques [7-21]. The geometry of 
the electrode makes it more difficult to analyse analytically, due to the discontinuity in the boundary 
conditions at the edge of the electrode. Oldham [22] accounted for the edge effects to obtain an 
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asymptotic expression for the short-time behaviour in the diffusion-limited regime. Subsequently, Aoki 
and co-workers developed analytical expressions using the Wiener-Hopf method for the short- [23] and 
long-time [24] chronoamperometric current responses in the same regime. Phillips and Jansons [25] also 
derived an analytical expression for the short-time current response using a Brownian motion model for 
the diffusion of the redox species, and, as part of a more general article on two-and three-dimensional 
diffusion in microelectrode chronoamperometry, Phillips [26] used matched asymptotic expansions to 
derive an expression for the long-time current response, which agreed with the first term in the series 
derived by Aoki et al. [24]. In both articles [25] and [26], the inlaid microband is a specific case of 
more general formulae that can also be applied to many other electrode geometries. Since the short-and 
long-time asymptotic expressions break down at intermediate times, Szabo et al. [27] and Aoki et al. [28] 
have developed approximate analytical expressions to bridge this gap. 

All of the expressions derived previously are only valid for the diffusion-limited regime, either due to 
extreme polarization of the electrode, or when the diffusion coefficients of the oxidant and the reductant 
are equal, D Q = Dr = D, and the following non-dimensional reaction rate (cf. [1]) is effectively infinite: 

kfd khd 

P=^- + -^- (2) 
Do D R 

To exploit the capacity of microelectrodes to access information about fast redox kinetics, it is invaluable 
to understand analytically the effect of the kinetics on the current response. It is also important that the 
theory accounts for unequal diffusion coefficients of the oxidant and reductant, since this inequality can 
be significant (see for example [29], where the diffusion coefficient ratio for hexaammineruthenium in 
different electrolytes is measured to be ~ 0.71). In this article, we include the effects of finite kinetics 
and unequal diffusion coefficients, and derive an asymptotic analytical expression for the long-time- 
dependent Faradaic chronoamperometric current response due to two-dimensional diffusion at an inlaid 
microband electrode. The solution is valid for time-scales such that d 2 / D <C t L 2 /D, where D is the 
typical size of the diffusion coefficients of the oxidant and the reductant, i.e., it is valid for time-scales 
much larger than the diffusive time-scale defined by the width of the electrode, but much shorter than 
the diffusive time-scale defined by the length of the electrode, when the assumption of two-dimensional 
diffusion ceases to be valid. For the reader who wishes to skip the mathematical derivation, the final 
expression for the long-time Faradaic current per unit axial length of the electrode (A-m^ 1 ) is detailed 
in Equation (57), and the details of the theoretical problem are displayed in Figure 2. The current 
depends on a function Q((3), which is found at each value of the parameter (3 (given by Equation (2)) 
by solving the integral equation detailed in Equation (41) and substituting the result into Equation (42), 
or equivalently by solving the boundary value problem given in Equation (70). This function must be 
found numerically, and we provide a working curve for Q(f3) for 10~ 2 < (3 < 10 3 in the Supplementary 
Information. We also derive explicit asymptotic formulae for Q{f3) when (3 is small or large. Comparison 
with the results of numerical calculations validates the analytical expression for the current response, and 
we demonstrate that the effect of unequal diffusion coefficients can be significant. Finally, we discuss 
the implications of using the Butler- Volmer framework [1] to model the rate constants of a one-step, 
one-electron redox reaction, and indicate that the width of the electrode must be chosen carefully to 
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access accurate information about the standard kinetic rate constant and the electron transfer coefficient 
from the long-time current response. 

2. Theory 

The theoretical problem is depicted in Figure 2. We consider two species, Ox and Red, diffusing 
in two-dimensions above an inlaid band electrode with its centre situated at x = 0, z = 0 and with 
width w — 2d (m). Throughout the article, tildes indicate dimensional entities. On the surface of 
the electrode, the species undergo the redox reaction described in Equation (1), and we assume that 
the potential of the electrode surface is held constant after an initial step, so that kf and kb are both 
constant. Neglecting any effects due to migration and convection, then the concentrations of Ox and Red, 
Co{x, z, t) and C R (x, z, t) (mol-m -3 ), satisfy the diffusion equation in z > 0 with unequal constant 
diffusion coefficients D Q and D R (m 2 -s _1 ). Hence 

DoV'Co^ 9 ^, D R V 2 C R = 9 ^, m~z>0 (3) 
ot at 

Initially the bulk concentrations of each species are constant everywhere: 

C 0 (x, z, 0) = C* 0 , C R (x, z, 0) = C* R , mz>0 (4) 

We assume that the bulk concentrations remain undisturbed far away from the electrode as the redox 
reaction at the electrode progresses, so that the far-field boundary conditions are: 

C 0 -> C* 0 , C R -> C R , as x 2 + z 2 -> oo (5) 

On the insulating surface, there is no flux of both species: 

dC 0 dC R 



dz dz 



0, on | x |> d, z = 0 (6) 



On the electrode surface, the boundary conditions are determined by the redox reaction in Equation (1) 
and conservation of matter: 



Do q~ — kjCo — k b C R , 



Do 9 ^ - -n R ac " 



> on | x |< d, z = 0 (7) 



oz oz 

The Faradaic current per unit axial length of the electrode, I(t) (A-m -1 ), is given by 

— -(x, 0, t)dx (8) 
d ® z 

where F is Faraday's constant (96,485 C-mol -1 ) and we recall that n is the number of electrons 
transferred in the redox reaction. 
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Figure 2. Schematic of the theoretical two-dimensional diffusion problem to be solved 
for the Faradaic current through an inlaid microband electrode. Two redox species with 
concentrations Co and Cr (mol-m~ 3 ) diffuse with unequal diffusion coefficients Do and 
D R (m 2 -s _1 ) above the electrode. Initially the bulk concentrations are uniform everywhere 
and equal to C* Q and C R ; we assume these remain undisturbed as x 2 + z 2 — > oo. On the 
surface of the electrode, | x \< d, z — 0, the boundary conditions model the redox reaction 
in Equation (1) and conservation of matter. The Faradaic current per unit axial length, I(t) 
(A-m -1 ), through the electrode is proportional to the integral of the diffusive flux at the 
electrode surface. 



x 



As x 2 + z 2 -» oo : C 0 -» C 0l C R -» C* R 



D 0 V 2 C 0 



dCp 
dt 



DrV 2 C r 



dC 



R 



dt 



x = —d 



x = d 



z = 0 



8C 0 dC R 
Insulating surface: — = ^ = 0 



x \> d, z = 0 



dz dz 



Electrode surface: Dq 



x \< d, z = 0 



dCp 
dz 



kfC 0 - k b C R 



dz 



dz 



dCo 

Faradaic current per unit axial length: lit) = —nFDg I — ^ — (x,0,t)dx 

-d oz 



To non-dimensionalise the problem, we choose the following scalings, where d is half the width of 
the band: 

7 ~ ~ d 2 

x = dx, z = dz, t = —t (9a) 

D 



D 



o 



DD 



o- 



D 



R 



DD, 



d 



D °l- 
-^-k b 

d 



(9b) 



a 



o 



C 0 - 



k f C 0 - k b C R 



Co, C_ 



R 



C R - 



k fCh - hC R 



c 



R 



(9c) 



nFD Q | kf ° h ~ hC * R I / 



(9d) 



Here D is the typical size of the diffusion coefficients of the oxidant and the reductant, Do and 
D R , and (3 is given by Equation (2). We remark that we have scaled time with the diffusive time-scale 
associated with the dimensions of the width of the electrode; over time-scales of this order, the diffusion 
can be regarded as two-dimensional. 
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Then, in terms of the non-dimensional variables, the problem becomes: 

BC BC 

D 0 V 2 C 0 = D R V 2 C R = mz>0,t>0 (10a) 
at Bt 

C 0 , C R = 0, in z > 0, t = 0 (10b) 

Co, C R -)> 0, as x 2 + z 2 -» oo (10c) 

¥^ = ^ = 0, on|*|>l,z = 0 (lOd) 



o 



Do dz ' on | x |< 1, z = 0 (lOe) 

kfCo — k h C R — /3, 



where (3 is given by Equation (2). The dimensionless Faradaic current per unit axial length through the 
electrode is given by 

/l ari 
-^(x,0,t)dx (11) 

2.1. Asymptotic Solution for the Long-time Transient Behaviour 

To solve this problem, we take the Laplace transform in time of the expressions detailed in 
Equation (10) to obtain: 

D 0 V 2 Co = sC 0 , D R W 2 C R = sC R , inz>0 (12a) 

Co, C R -»• 0, as x 2 + z 2 -»■ oo (12b) 

BC 0 BC R i ■ n . 10 v 

— — = — — =0, on \ x \> 1, z = 0 (12c) 



_^=J D 0 Bz ' on | .t |< 1, z = 0 (12d) 

l^/c/Co — k b CR — — , 

The long-time behaviour, t 1, can be deduced from the small s behaviour, |s| <C 1, in Laplace 
transform space; in dimensional terms, this corresponds to time-scales such that t ^> d 2 / D. Since the 
term f3/s in the redox reaction boundary condition on the electrode, detailed in Equation (12d), blows 
up as s 0, we re-scale the problem by setting Co,r = co,r/s, which gives: 

D 0 V 2 c 0 = sc 0 , D r V 2 cr = sc R , in z > 0 (13a) 

Co, cr — > 0, as x 2 + z 2 — > oo (13b) 

dc 0 dc R II -, n n ~ , 

— — = — — = 0, on \ x \> 1, z = 0 (13c) 

Bz Bz 

DrBcr 



Bcq 
Bz 



D 0 Bz' on | x |< 1, 2 > 0 (13d) 

kfCo- k b c R - (3, 



We use the method of matched asymptotic expansions to solve the problem given by Equation (13) 
approximately for | s |-C L In an inner region near the electrode, over 0(1) length-scales, the O(s) 
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terms on the right-hand side of the governing equations in Equation (13a) can be neglected. However 
these terms cannot be neglected in an outer region when the length-scales are 0(s -1 / 2 ), such that 

X = X/yfe, z = Zj4~s (14) 

In this outer region, the terms on both sides of the governing equations in Equation (13a) are of 
similar size. The method of matched asymptotic expansions requires us to search for solutions in the 
inner and outer regions, and then to match them in an intermediate region where the two solutions have 
an overlapping region of validity. 

In the inner region, we denote the approximate solutions by ho, R (x, z, s), and they satisfy Laplace's 
equation and the boundary conditions on the surface of the electrode and the insulator: 



S/ 2 h 0 = 0, V z h R = 0, in z > 0 (15a) 

dhn dhn 

—^ = —^ = 0, an he > 1,* = 0 ( 15b ) 

oz oz 

D R dh R 



dh 



o 



dz 



D 0 dz ' on | x |< 1, z = 0 (15c) 

kfho - k b h R - (3, 



In the outer region, we denote the approximate solutions by Ho, R (X, Z, s), and they satisfy the 
modified Helmholtz equations: 

D 0 V 2 H 0 = Ho, D R V 2 H R = H R , in Z > 0 (16a) 

with boundary conditions 

H 0 , Hr 0, as X 2 + Z 2 ->■ oo (16b) 
Then the problem is solved approximately by matching the inner solutions, ho, R , to the outer solutions 

Ho,r- 

The solutions to the inner problem given by Equation (15) have the following form: 

_ 0 (x z) 

h 0 (x,z,s) = A 0 + J2r^rA (17a) 



h R (x, z,s) = — 
l> r 



(17b) 



where A 0 and B 0 are constants, and the solutions for and Ri, that allow the conservation of matter 
boundary condition in Equation (15c) to be satisfied, can be written using a Green's function as: 

1 f 1 dO, 



Oi{x, z) -- 

Z7T 



-^(u,0)\og((x-u) 2 + z 2 ^du + A i (18) 
Ri(x,z) = -—J -^(u, 0)\og((x-u) 2 + z 2 ^)du + B l (19) 



where again A t and Bi are constants. The redox reaction boundary condition on the electrode surface 
given by Equation (15c) implies that: 

k f A 0 - k b ^B 0 = 0 (20) 



Sensors 2013, 13 



633 



and that the following integral equations must be satisfied on the surface of the electrode for i > 1: 

dO Q F dO D 

—i(x. 0) = - / — Mu, 0)log I x-u I du + kfAi - k b -^-Bi, for I x |< 1 (21) 
02 7T J_ x dz D R 

The solutions to the outer problem in Equation (16) are given by (writing R 2 = X 2 + Z 2 ): 

Ho{X,Z,s) = j2<*K 0 (*) (22a) 



^io g (vVr \VD 0 
h r (x,z,s) = ^itv^rv^i^r) (22b) 

where Cj and Dj are constants, and K 0 (-) is a modified Bessel function of the second kind [30]. We 
remark that the branch cut of y/s in the relation expressed in Equation (14) must be taken along the 
negative real axis, so that the vanishing boundary conditions given by Equation (16b) at infinity are 
satisfied. The unknown constants in the inner and outer solutions, Aj, Bi, Cj and D L , must be found by 
asymptotic matching. 

We match the inner and outer solutions using the intermediate variable technique [31]: 

n = r^s a = RV~s~ {l - a) (23) 

where r 2 = x 2 + z 2 and 0 < a < 1. Written in terms of the intermediate variable 77, the functions 0« and 
Ri in the inner solution detailed in Equation (17), and given by the expressions in Equations (18) and 
(19), become: 

1 „ . ^ f 1 dO, 



1 f dO 

Oi = - (log r/- a log v^) / -^(u, 0) du + A { + OW?) (24) 

71 J_ 1 OZ 

X f ^ dO 

Ri = (log 77- a log v^) / — 0) dw + 5, + 0(v^ Q ) (25) 

whilst the modified Bessel functions in the outer solutions expressed in Equation (22) become: 

K 0 (-^=A = - 7 -(l-a)logv^-logr7 + log v / ^ + log2 + 0( S (1 - cl) log S ) (26) 

K 0 (-^=\ = -7- (1 -a)logv^-logr7 + log V / ^ + log2 + 0(s (1 - Q) logs) (27) 

where 7 is Euler's constant. 

In the intermediate region, therefore, the inner expansion can be written: 

- A a f 1 dO Xt 

h 0 ~ A 0 / — — (u, 0)du 

7T y_! dz 

- DM ocD 0 f 1 dO x 

h R — + — — / - s — (u, 0 dw 

Dq\^ ( a f 1 dO i+ i 



£ ("~ £ («, 0) d« + ^ £ 0) d« - 2*) (log (28b) 
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and the outer expansion becomes: 



H 0 ~ -Ci(l - a) + £; (a(log2 - 7 - log?? + log v 7 ^) - (1 - a)C i+ i) (log 



(29a) 



-n r ~ — ( 1 — a 



oo 

(A (log 2 - 7 - logr/+log v / ^)-(l-«)A + i) (log (29b) 

j=i 



Comparing similar powers of (log y/s) 1 in Equations (28) and (29) gives at 0(1) 

a f 1 dOi 



A 0 - 
B 0 + 



7T J_ 1 8z 

a f 1 dO x 



7T J_ t 3Z 



(u, 0) du 
(u, 0) du 



-d(l-a) 
-D x {l-a) 



which implies that: 



An = -fin 



-Ci = A 



i /- 1 ad 



7T J_i 8z 



u, 0) dw 



Using the condition given by Equation (20), we find that: 

1 f 1 d0 1 



7T J_ 1 8z 



u, 0) du = 1 



so that 



A 0 = -B 0 = -Ci = A = 1 
At 0(log v^) * (* — 1,2,3,.. .), the matching process between Equations (28) and (29) gives: 



(30) 
(31) 

(32) 

(33) 
(34) 



a f 1 dO i+1/ log// f 1 dOi 

u, 0) du- 1 



7T 



(_1 9-2 

a f 1 dO i+ i 



[u, 0) du- 



ix J_ x dz 
log 7? f 1 dOi 



7T J_ 1 dz 7T J_ 1 dz 

Equations (35) and (36) imply firstly that: 



[u, 0) du + A i = d (log 2-7 - log r?+log y/D 0 )-( 1-a) C i+ i 

(35) 

u, 0) = A(log 2-7 - log 77+log v / ^)-(l - «) A+i 

(36) 



d = -A = — 



1 f 1 dd 



7T 



'it, 0) du, for % > 1 



(37) 



and, secondly, that: 



C m = C i (log2-7 + log v / A^)-A, 
A+i = A (log 2 - 7 + log a/Ar) - Bj, 



for i > 1 



(38) 



Now we note that the conditions given by Equations (32) and (37) provide solvability conditions for the 
integral equations in Equation (21), namely that: 

1 f 1 dd 



7T J _ x dz 



(u, 0) du = —C i} for i > 1 



(39) 



Sensors 2013, 13 



635 



Writing 



^ v x, 0) = /T 1 ^ ; A t - fc^ij ^(z) (40) 
then the integral equations in Equation (21) and solvability conditions given by Equation (39) become: 

F(x) = - [ F(u) \og\x-u\ du + (3, for | x |< 1 (41) 



7T 

with 



-i 



where 



l£/ {u)iu= w) (42) 

Q(/3) = -^''C- 1 (k f A, - h^B,] (43) 

Solution of the integral equation given in Equation (41) for F(x) and substitution into Equation (42) 
determines Q{0) as a function of /3, and then Equation (43) gives another relationship between Ai, Bi 
and Ci. 

Eliminating Ai, B^ and Di from Equations (37), (38) and (43), we find a recurrence relation for C;: 



C m = |log2-7 + g(/3) + log v / ^+^r 1 ^logy / ^ja (44) 
with initial condition C\ — — 1 from Equation (34). This has solution for z > 1 given by: 

Ci = - (k)g 2 - 7 + Q{fi) + log + fc^- 1 ^ log (45) 

Finally, the Laplace transform of the current per unit axial length, I{s), is given by the following: 

I(s) ~ -y ^(«,0, s)du (46) 



oo 

7T 



^a(logVi)^ (48) 



s 

where we have used the expression for h Q given in Equation (17a), and have employed the relationship 
given in Equation (39). Substitution of the solution for Cj detailed in Equation (45) into Equation (48) 
gives: 

71 V" ( i_o i ^vm , i_ /FT" i i. , [D~ 



7 ( s ) ~ -^^log2-7 + Q(/3) + logv / ^ + ^- 1 ^log^/ 7 ^) (logVi)- (49) 

"7T\ 1 



log ^ - (log 2 - 7 + W) + log VD^ + k b ^ log ^ 



(50) 
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This can be re- written in similar form to the solution in Phillips [26]: 



where <p is given by: 



/W ~-'t) ,og( S -.) + Ulog Do (51) 

4> = 2 (log 2 - 7 + Q(/3)) + log (52) 
Then Phillips [26] shows that this can be inverted to give the long-time current as: 

I(t) ~ -lixT^Dot) (53) 

where 

1(a) = r , x du (54) 



o u^(logu) 2 + 7T 2 j 



Integration by parts (to remove the integrable singularity at u = 0), and the change of variable v = au, 
casts Equation (54) into a more amenable form for numerical evaluation: 



1 1 f°° 

1(a) — - H — / e ^arctan 

2 ^ Jo 



-\og(v/a) 

7T 



dv (55) 



Phillips [26] tabulated values of the function 1(a) and showed that it can be expanded asymptotically in 
inverse powers of the logarithm for large a, such that 

1(a) = (log a)" 1 - 7 (loga)- 2 " " ^ {^ga)^ + OQoga)" 4 (56) 

where 7 is Euler's constant. Although this expression is simpler to implement, Phillips [26] emphasises 
that the expression in Equation (54), or equivalently Equation (55), is more accurate as it incorporates 
all the logarithmic terms, so that the error in the current is algebraic rather than logarithmic. 

In dimensional terms, the expression in Equation (53) for the Faradaic current per unit axial length 
(A-m -1 ) through the microband electrode becomes: 

fc* 0 -hc* R \ 

I(t) ~ -2imFD 0 I j l^Dot/d 2 ) (57a) 



k f D 



f^R 




where !(•) is evaluated from Equation (55) and 0 is given by: 

0 = 2(log2- 7 + Q(/3))+ (^1 + 2^ ] log(^) (57b) 

We recall that (3 is given by Equation (2), and the function Q(f3) in Equation (57b) must be calculated 
numerically by solving the integral equation in Equation (41) for F(x) and substituting the result into the 
expression in Equation (42). In order that the experimentalist does not have to perform these numerical 
calculations themselves, we have provided a working curve for Q(f3) over the range < (3 < 10 3 as 
the file "QJoeta_working_curve . csv" in the Supplementary Information. The calculation of this 
working curve is detailed in the next section, where we also provide asymptotic analytical expressions 
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for small and large 8, given by Equations (62) and (63), which can be used to calculate Q(/3) outside the 
range of the working curve. 

We remark that the Faradaic current per unit axial length through a laminar electrode in free space is 
simply double the expression in Equation (57a), since the current is identical through both sides of the 
lamina. 



2.2. Calculation of Q {(3) 

It is useful for the experimentalist to have a working curve to evaluate the function Q(/3) in Equation 
(57), and in this section we describe how we have calculated it. We also provide analytical asymptotic 
expressions that can be used to calculate Q{8) for small and large 8. 

For a particular value of 8, Q(8) is calculated by solving the integral equation in Equation (41) for 
F(x) and inserting the result into the expression in Equation (42). We solved the integral equation by 
removing the logarithmic singularity and applying the Nystrom method (cf. Delves and Mohamed [32]). 
Consider N Gauss-Legendre quadrature weights, W{, at the abscissae X{ in the interval (—1, 1). At each 
point X{, the logarithmic singularity in the integral equation detailed in Equation (41) can be removed by 
writing it as: 



F (?i) = ^J i F ( s ) ~ F ( x *)) log | a* - s | ds 



+ ^-((l-x i )\og(l-x i ) + (l + x i )log(l + x i )-2)F(x i )+8, for % = 1, . . . ,N (58) 

7T 

The integrals in Equation (58) can be evaluated using Gauss-Legendre quadrature at the same abscissae 
to give: 

8 l ~ l (3 N 

F(x i ) = - S ^w k (F(x k ) - F(xi)) log | Xi-x k | H — ^ w k (F(x k )-F(xi)) log | x^ - x k \ 

k=l k=i+l 

+ -((1 - log(l - Xl ) + (1 + Xl ) log(l + Xi) - 2)F(xi) + (3 (59) 

7T 

with appropriate care taken whenever % = 1 or i = N. This represents N linear equations to be solved 
for the iV unknowns F(xi). Then Q(f3) can be calculated from Equation (42) as follows: 

For each B, we calculated Q(B) using = 400 Gauss-Legendre weights. The values of 8 used to 
calculate the points on the working curve for Q(8) were logarithmically spaced between 0.01 and 1000 
as follows: 

0; = io-^/ 2000 , j = 0, . . . , 2000 (61) 

Convergence was slower for larger 8j, and the maximum absolute change in Q((3j) on doubling iV from 
200 to 400 was less than 10 -6 at 8j = 1000. The working curve is attached in the Supplementary 
Information as a CSV file named "Q_beta_working_curve . csv". 
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Asymptotic expressions can also be derived for small and large (3. For small (3, the function Q((3) 
takes the following form: 

Q(/3) = ^ + ^-log2 + 0(/3) (62) 
For large f3, the asymptotic solution is given by (proof in the appendix): 

Q(/3) = log 2 + \ log/3 + ot/T 1 log/3) (63) 
ix p 

We include a log-log plot of Q{(3) in Figure 3, which shows that the agreement between the asymptotic 
solutions (dashed lines) and the numerically calculated values (solid lines) is good. For (3 < 0.01, the 
expression in Equation (62) can be used to calculate Q((3), while for (3 > 1000, the expression in 
Equation (63) provides a reasonable approximation. For all other values, the numerically calculated 
working curve can be used. 

Figure 3. Log-log plot of the function Q{(3). The solid line has been calculated numerically 
as described in Section 2.2, and the dashed lines have been plotted using the asymptotic 
expressions given by Equation (62) for (3 <C 1, and Equation (63) for (3 ^> 1. 
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3. Results and Discussion 

In the previous section, we have derived the asymptotic solution detailed in Equation (57) for the 
long-time-dependent chronoamperometric current per unit axial length due to two-dimensional diffusion 
at an inlaid microband electrode. The solution allows for finite kinetics and unequal diffusion coefficients 
of the oxidant and the reductant. In this section, we show that the expression in Equation (57) reduces 
to the first term in the series derived by Aoki et al. [24] (see also Phillips [26]) for the current in 
the diffusion-limited regime due to extreme polarization. We also present the simplified form of 
Equation (57) for reversible reactions; in the case of identical diffusion coefficients, this expression 
also reduces to the first term in the series derived by Aoki et al. [24]. We compare the analytical 
solution to the results of numerical calculations, and we illustrate that unequal diffusion coefficients 
can cause significant differences in the current response. Finally, we discuss the formula for the current 
response due to a one-step, one-electron redox reaction whenever the rate constants are modelled by the 
Butler- Volmer expressions [1], and we indicate that the width of the electrode must be chosen carefully 
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to allow accurate estimates of the standard kinetic rate constant and the electron transfer coefficient to 
be obtained from the long-time current response. 

3.1. Diffusion-limited Currents due to Extreme Polarization 

The diffusion-limited reduction current per unit axial length due to extreme polarization corresponds 
to kf — > oo and k b — > 0. This means that (3 — > oo, and the expression in Equation (57) for the current 
response becomes: 

lit) ~ -2irnFD 0 C* 0 iflQe-^Dot/d 2 ^ (64) 

where we have used Equation (63) to set Q(oo) = log 2. This expression agrees with first term in the 
series derived by Aoki et al. [24] (see also Phillips [26]). 

Similarly, the diffusion-limited oxidation current per unit axial length due to extreme polarization 
corresponds to kf — > 0 and k b — > oo, which gives: 

7(f) ~ 2nnFD R C* R l(lQe- 2 ~<D R t/d 2 ^ (65) 

As expected, this is equivalent to Equation (64) with the sign changed, and the subscripts O and R 
interchanged. 

3.2. Reversible Reactions 

For reversible reactions, kf, k b — > oo such that kf/k b = 0(1). In this case, since /? — > oo and 
Q(oo) = log 2, we obtain: 

(C* 0 - ^C* R \ N 
I(t) ~ -2-imFDo l(e^D 0 t/d 2 ) (66a) 

\ 1 + k^ D^ J 

where 

In the case of identical diffusion coefficients, this expression also reduces to the first term in the series 
derived by Aoki et al. [24], since the logarithmic term in Equation (66b) vanishes. 

3.3. Comparison with Numerical Calculations 

To calculate the current numerically, we employed the fully implicit finite-difference method (FIFD) 
devised by Gavaghan [33,34] for the inlaid disk electrode, but adapted for the Cartesian coordinate 
system appropriate to the microband electrode. We simulated the non-dimensional problem detailed 
in Equation (10) on the finite domain 0 < x < x max = 201, 0 < z < z max = 200 for 0 < t < 
t max = 100. Considerations of symmetry allowed us to simulate only half the domain (x > 0) with 
extra symmetry boundary conditions dCo/dx = dCn/dx = 0 on x = 0. The finite boundaries of 
the domain, x max = 201 and z max = 200, were chosen to satisfy the criterion that x max , z max > 
max{6\/ D 0 t mSLX , 6v 'D nt mSLX } (cf. [17]); this ensures that application of the conditions C Q = Cr = 0 at 
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the finite boundaries x max and ,2 max , instead of infinity, does not distort the current through the electrode. 
Gavaghan's method uses a spatial grid that expands exponentially away from the edge of the electrode 
in order to resolve accurately the large flux in the neighbourhood of the edge. Here we used the same 
grid parameters as suggested in [33,34]: the spatial steps in the x- and ^-directions nearest the edge of 
the electrode at x — 1 are given by h last = 8 x 1CT 5 , and the expansion factor is given by / = 1.175. 
The initial time-step was taken to be 1CT 5 , and was increased by a factor of 10 every thousand steps. 

Figure 4. Plots of I(t)//3 versus t for different rate constants: (a) kf = 1, k b = 1; 
(b) kf — 5, kf, — 1; and (c) kf = 1, k b = 5. The solid lines are calculated using the 
approximate analytical expression in Equation (53) for the long-time current response, while 
the symbols are plotted from the results of the numerical calculations described in Section 
3.3. The symbols correspond to different ratios of the diffusion coefficients: Do = 1, 
D R = 0.5 (triangles), D 0 = 1, D R = 1 (squares), and, D 0 = 1, D R = 2 (circles). It 
is apparent from the plots that the effect of unequal diffusion coefficients can be significant. 




0 20 40 60 80 100 



t 

(c) kf = 1, kb = 5 



We plot I(t) / (5 in Figure 4 (a-c) for 0 < t < 100 for different rate constants: (a) kf = 1, k b = 1, (b) 
kf = 5, k b = 1, and (c) kf = 1, k b = 5. We have divided the current by (3 to capture completely the 
effect of varying the ratio of the diffusion coefficients (cf. the factor of in the non-dimensionalisation 
of the current in Equation (9d)). We note that (3 (given by Equation (2)) can be written in terms of the 
non-dimensional variables as: 

f3 = kf + h 1 ^. (67) 

The different symbols in each plot correspond to different ratios of the diffusion coefficients: Do = 1, 
D R = 0.5 (triangles); D Q = 1, D R = 1 (squares); and, D Q = 1, D R = 2 (circles). The solid 
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lines are calculated using the analytical expression in Equation (53) for the non-dimensional long-time 
current response, with 0 given by Equation (52) and the integral 1(a) evaluated using the expression in 
Equation (55), while the symbols are plotted from the numerical results. It is apparent from 
Figure 4 (a-c) that unequal diffusion coefficients can substantially alter the magnitude of the current 
response. In Figure 5 (a-c), we plot the percentage difference between the numerical results and the 
analytical expression as a function of time for the same parameters. As expected the analytical expression 
approaches the numerical results more closely for longer times, since this is when it becomes valid. For 
t > 5, the maximum difference is less than 1.75% for all the parameters considered. (We also remark 
that comparison with the reference values calculated by Britz et al, [18] for the diffusion-limited current 
indicates a similar level of accuracy; for t > 5, the percentage difference in this regime is less than 1%.) 

Figure 5. Log-linear plots of the time- varying percentage difference between the analytical 
expression given by Equation (53) for the current response and numerical results calculated 
as described in Section 3.3, using the same parameters as the plots in Figure 4: (a) kf = 1, 
kb = 1, (b) kf — 5, kb — 1, and (c) kf = 1, kb = 5. The symbols correspond to different 
ratios of the diffusion coefficients: D Q = 1, D R = 0.5 (triangles), D Q = 1, D R = 1 
(squares), and, Do = 1, Dr = 2 (circles). The analytical expression approaches the 
numerical results for longer times, when it becomes valid. For t > 5, the maximum 
difference is less than 1.75% for all the parameters considered. 




(c) kf = 1, kb = 5 

3.4. Long-time Current Response for a One-step, One-electron Redox Reaction Modelled in the Butler- 
Volmer Framework 

So far we have not assumed any model for the rate constants of the redox reaction, kf and k b . For 
a one-step, one-electron process, so that n — 1 for the redox reaction detailed in Equation (1), the 
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Butler- Volmer model [1] describes the rate constants as an exponential function of the potential E (V) 
at the electrode surface: 

k f = ~k oe - a w^°\ k b = ~k 0 e {1 - a) w(*-^ (68) 

where k 0 (m-s -1 ) is the standard kinetic rate constant, 0 < a < 1 is the electron transfer coefficient, and 
E 0 (V) is the formal oxidation potential. The remaining parameters are Faraday's constant, F (96,485 
C-mol -1 ), the universal gas constant, R (8.3145 J-K _1 -mol _1 ), and the temperature, T (K). Substitution 
of this form of kf and kb into Equation (57), and setting n — 1, gives the Butler- Volmer long-time 
current response as: 

I(t) ~ -IvFDo I ° * I Z( e^Dot/d' ) (69a) 



where 



2 (log 2 - 7 + W)) + ( 1 + e-w^-^S^ ) log ( ^ ) (69b) 




and (3, given by Equation (2), becomes: 

-aX(B-Bo) , Al-a)-L{E-E 0 ) 





It is interesting to observe from Equation (69) that the only impact of the standard kinetic rate constant, 
ko, and the electron transfer coefficient, a, on the long-time current response occurs through the function 
Q((3). The sensitivity of Q{(3) to the parameters k 0 and a depends on its gradient, which tends to zero 
as (3 — > oo (see the asymptotic expression in Equation (63) and Figure 3). In addition, as f3 — > 0, 
the current response, and its sensitivity to ko and a, also tends to zero. This suggests that there is a 
range of (3 over which the sensitivity to the parameters k 0 and a is significant enough to allow them 
to be estimated accurately using the long-time current response; outside this range the sensitivity to 
the parameters will be too small. Since the size of f3 can be controlled by changing the width of the 
microband d (cf. Equation (69c)), this implies that the size of the electrode should be chosen such that (3 
lies within this sensitive range. However, as this range will depend on all of the system parameters, the 
appropriate size of d can only be estimated before the experiment using a priori estimates of each of the 
parameters. It can then be verified a posteriori that the sensitivity of the current response was significant 
enough to provide accurate estimates of the parameters, e.g., by calculating confidence limits. We intend 
to explore the choice of electrode size more carefully in a future article. 



4. Conclusions 



We have derived the approximate asymptotic expression given in Equation (57) for the long-time 
chronoamperometric current response at an inlaid microband electrode whose axial length is much larger 
than its width so that the diffusion is effectively two-dimensional. For a laminar electrode in free space, 
the current response is simply double the expression given in Equation (57). The solution allows for finite 
redox kinetics and unequal diffusion coefficients of the oxidant and reductant. As an input, it requires 
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calculation of a function Q(/3), which we have determined numerically and provided as a working curve 
in the Supplementary Information. We have also provided simple asymptotic expressions that can be 
used to calculate Q(/3) whenever (5 is small or large. The expression for the current response has been 
validated by comparison with the results of numerical calculations, and we have demonstrated that the 
effect of unequal diffusion coefficients on the current response can be substantial. Finally, we have 
discussed the form of the long-time current response due to a one-step, one-electron redox reaction 
when the rate constants are modelled using the Butler- Volmer framework, so that the current response 
is given by the expression in Equation (69). The formula indicates the importance of choosing the width 
of the microband carefully to ensure accurate estimation of the standard kinetic rate constant and the 
electron transfer coefficient from the long-time current response. We intend to investigate this more 
thoroughly in a future article. 
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Appendix 

A. Large-/? Asymptotic Solution for Q{0) 

To derive the large-/? asymptotic solution for Q(f3), we employ a similar argument to that used by 
Phillips [35] to determine the asymptotic solution for near-diffusion-limited currents at inlaid electrodes 
when the diffusion is three-dimensional. 

As described in the main body of the article, the function Q(f3) can be found by solving the integral 
equation detailed in Equation (41) and substituting the result into Equation (42). However, to proceed 
with the large-/? asymptotics, we consider an alternative, but equivalent, formulation of the problem. The 
function Q{(3) can be also be determined by finding the constant coefficient as r — > oo (which is equal 
to Q{0)) in the solution to the following boundary value problem: 

in z > 0 (70a) 

on | x |> 1, z = 0 (70b) 

on | x |< 1, z = 0 (70c) 

as r — > oo (70d) 

As (3 — > oo, the leading-order solution to the problem in Equation (70), denoted by Goo, satisfies: 

inz>0 (71a) 

on | x |> 1, z = 0 (71b) 

on | x \< 1, z = 0 (71c) 
as r — > oo (7 Id) 

with solution derived using conformal mapping (see [26]) to be: 

Goo = Re\og(x + Iz + ((x + Iz) 2 -1)^ (72) 

Finding the constant coefficient in this solution as r — > oo implies that: 

Q(oo) = log2 (73) 



V 2 G 

dG 

dz 
dG 

dz 



0. 

o, 



G ~ log r, 



V'Goo = 0, 
dG^ =Q 

dz 

Goo — 0, 



Goo ~ log r, 
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To find higher order corrections to this expression, we employ a similar technique to Phillips [35], and 
use Green's second identity to write: 

j\l-G)^-(l-G OB )^dx = (l-G^-il-G^dl (74) 

where G is the solution to the entire problem given in Equation (70). Here the integral on the left-hand 
side is along the surface of the electrode on z = 0, while, on the right-hand side, the integral is around 
a semi-circle of radius R in z > 0 centred on the origin (denoted by T R ) and d/dn denotes the outward 
normal derivative. Now, the condition that G, (?„ ~ log r as r — > oo implies that: 

ax = 7T, / — — dr = 7r (75) 



l dz J_i dz 

Also, on the surface of the electrode: 

1 dG 

Goo = 0, G = - — , on|x|<l, z = 0 (76) 

p dz 

and, in the far-field: 

Goo ~ log r + log 2, G ~ logr + Q{/3), as r ^ oo (77) 
Substituting the expressions in Equations (75)-(77) into Equation (74), and letting R — > oo, we find that: 

Q( /3 )=lo g2+ i 3 /_ 1 J^d, (78) 

The leading-order solution for G is G^, and this breaks down as a singular perturbation in an 0(/3 _1 ) 
region of x — ±1. Taking 1 ^> 5 ^> and using symmetry about x = 0, we can write the integral on 
the right-hand side of Equation (78) as: 

f 1 9G 9G - d. ~ 2 f ' 2 dx + 2 T ^ d* (79) 

oz oz 1 i.i. 

Since 



j <9z <9z 7 0 \ dz J Ji_ s dz dz 



dG 



oo 



:, on \x\<l, z = 0 (80) 



•92 Vl - 

the first integral on the right-hand side of Equation (79) can be evaluated to be: 

2 dx = \og2-\ogd + 0(5) (81) 



J o V / 

Substituting x = \ — (3~ l u into the second integral on the right-hand side of Equation (79), we obtain: 

2 r 1 dGdG^ dx ^^ f SPd ^J_ du (g2) 

where 5(5 3> 1. But as u — Y oo, G ~ Goo, so that for u > M for some large constant M: 

<9G dG^ yft 
•72; dz J2u 
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Therefore 

2 f SfS 9G 1 J [ sp 1 , 2 /- M 9G 1 , 

dw ~ / - du + -= / — = dit (84) 



V/3 io Jm u \ffi Jo 9z y/2u 

The contribution from the first integral on the right-hand side dominates, so that: 

dGdG r 



lis dz dz 

Adding the two contributions from Equations (81) and (85), we find that: 



dx ~ log 5 + log ^ (85) 



nf3 J_i dz dz 7l(3 
and hence, from Equation (78), that: 



1 f 1 dGdG^ 1 

dx ~ — logp, (86) 



Q(/3)~log2+^-log/3 (87) 

7Tp 
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